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Відновлення потенціалу сили тяжіння за значеннями модуля його градієнта в задачі Алексідзе 


Нелінійна гранична задача Алексідзе для рівняння Лапласа постулює аналітичне продовження сили тяжіння 
в глобальних областях. Її розв'язок як потенціал простого шару на поверхні Ляпунова є послідовністю розв'яз- 
ків зовнішніх граничних задач Неймана для рівняння Лапласа, якщо розв'язок не дуже ухиляється від заданого. 
Густина простого шару визначається з інтегрального рівняння Фредгольма. Вказані умови коректності. 
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Актуальність проблеми. Вирішення прикладних завдань гравіметрії і геодезії, пов'язаних з вивченням фі- 
гури, внутрішньої будови Землі, чи проявів її зовнішнього гравітаційного поля, потребує відомостей про розпо- 
діл значень потенціалу сили тяжіння чи модуля його градієнта. Для тлумачення даних аномалій потенціальних 
полів мало створити надійну методику побудови на їх основі геологічно змістовних моделей глибинної будови 
земної кори; вона має бути адекватною вимогам сьогодення. Методи видобутку з накопичених геофізичних 
даних всієї повноти необхідної інформації недостатньо розвинуті. Чинні методи трансформації потенціальних 
полів Г1, 2) завершують теорію потенціалу і інформативні у вивченні локальних особливостей приповерхневої 
будови планети: вони дають змогу розв'язувати з достатньою точністю ту чи іншу задачу лише за умови за- 
дання вхідної інформації в локальних областях певної малої міри. Спроби стикувати результати локальних 
розв'язків при відновленні інформації в глобальному масштабі зазнають невдачі |З) через відсутність точних 
граничних даних для розв'язання відповідних граничних задач - Діріхле, Неймана, Стокса-Молоденського - 
для рівняння Лапласа. Немає змоги і прямо вимірювати значення гравітаційного потенціалу, зате доступні да- 
ні гравімагнітних спостережень, що є значеннями приростів модуля градієнта потенціалу сили тяжіння (МГП- 
СТ). Варто скористатися ними при розробці схем трансформацій потенціалу в глобальній області. 

Серед трансформацій потенціальних полів чільне місце посідає задача наближеного аналітичного продов- 
ження потенціалу сили тяжіння: в моделі геологічного середовища наближено задано або диференціальний 
оператор, або граничні умови. Найчастіше використовують модель з точним диференціальним оператором 


(гармонічна апроксимація значень сили тяжіння 2, (Х) в задачі Молоденського |41. Відновити з гарантова- 


ною точністю поле сили тяжіння у зовнішньому просторі за його гармонічним наближенням можна шляхом 
розв'язання відповідної граничної задачі в області малої міри, оскільки у вищезгаданих граничних задачах да- 
ні спостережень слугують лише наближеними граничними умовами через негармонічність відповідних функ- 
цій. Через негармонічність оператора трансформанти |5) (принаймні, на суходолі) різко зростає похибка виз- 


начення потенціалу з, (х) зі збільшенням розміру локальної області. Трансформації на основі цих задач ма- 


ють гарантовану точність лише в областях малої міри (як правило, 19х19 і з точністю до невизначеної сталої, 
залежної від геометрії області). Гарантована точність розв'язку суттєво залежить від геометрії (рельєфу і роз- 
мірів) локальної області |6Ї, а критерії для поєднання локальних розв'язків в глобальних побудовах відсутні. 





Рис. 1. Розходження векторів реальної і нормальної сили тяжіння. 
У класичному способі визначення гравіаномалій 2. (У) не враховують те, що в точках земної поверхні ве- 


ктори реальної і нормальної сили тяжіння можуть мати різні напрями через розходження поверхні земного ре- 
льєфу та референц-еліпсоїда (рис. 1), що особливо яскраво проявляється на стику континент-океан і гірських 
районах. Просторову орієнтацію гравіметрів (нахили приладів, зумовлені різним ступенем кривизни еквіпотен- 
ціальних поверхонь, що проходять через задані пункти вимірювань) характеризує приріст кута 


ад, со5(п, , т,)- со5(п,. й та) між нормалями до земної та еквіпотенціальної поверхонь в точці вимі- 
рювань. Ні цю, ані жодну іншу величину, скажімо, значення напрямних косинусів СО5(п, Х, ), і -- 1,3, що виз- 
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начають напрям сили тяжіння, не вимірюють через складність організації таких спостережень в польових умо- 
вах. Розходження (приріст) цих напрямів в реальних умовах коливається від 0" до 40" і може призвести до то- 
го, що глибинні неоднорідності, які зумовлюють ці розходження, не відображатимуться в класичних аномаліях. 

Способи розв'язання. Необхідність вивчення нелінійної граничної задачі відновлення потенціалу сили тя- 
жіння за значеннями МГПСТ продиктована практичною непридатністю класичних схем відновлення потен- 
ціалу в глобальній області. При переході до глобальних побудов густинних моделей земної кори на основі да- 
них регіональних спостережень слід розв'язати задачу аналітичного продовження аномалій сили тяжіння, 
сформульованої в (ЗІ. В праці |6) для цього окреслено дві альтернативи, виходячи з аналізу характеристичних 
властивостей МГПСТ. 

Перша з них - відшукання такого диференціального оператора, що анулює значення МГПСТ поза областю 
розташування тяжіючих мас, та розв'язання для нього відповідної лінійної граничної задачі, зосібна, зовніш- 
ньої задачі Діріхле для лінійного диференціального рівняння типу Клейна-Гордона (яке моделює значення 
МГПСТ в області, не зайнятій тяжіючими масами) з невідомим змінним коефіцієнтом, що відповідає обраному 
нормальному потенціалу. Подібному рівнянню задовольняють і аномалії сили тяжіння |5). Цей спосіб ефектив- 
ний для продовження аномалій сили тяжіння, а у випадку продовження значень повного градієнта потенціалу 
не дає бажаних наслідків. В |З) обгрунтована і реалізована подібна схема у вигляді послідовності розв'язків 
задачі Неймана для рівняння Лапласа, які визначають збурювальний потенціал. 

Друга альтернатива -- постановка такої нелінійної граничної задачі для рівняння Лапласа, в крайових умо- 
вах якої безпосередньо задіяні значення сили тяжіння. Таку задачу вперше сформулював Алексідзе в роботі 
ІЗ1, а в І61 її переформульовано з граничними даними для класу поверхонь Ляпунова за умови, що відновлю- 
ваний потенціал не дуже відхиляється від заданого. 

Ради пошуку точніших способів відновлення наближень сили тяжіння довелось вийти за рамки задачі Ді- 
ріхле і шукати уточнення коефіцієнта рівняння, первісно обчисленого для нормального потенціалу, що при- 
звело до побудови послідовних наближень потенціалу за граничними значеннями МГПСТ. Перехід до задачі 
відновлення потенціалу усунув необхідність обчислення наступних наближень як коефіцієнтів рівняння сили 
тяжіння, так і самих значень сили тяжіння, оскільки останні тепер можна знайти не лише з розв'язання задачі 
Діріхле для рівняння сили тяжіння, а і (що простіше) з безпосереднього диференціювання відновленого по- 
тенціалу. 

У зв'язку з цим некласична задача гравіметрії про відновлення потенціалу за значеннями МГПСТ набуває 
особливої ваги в колі обернених задач теорії потенціалу. Один з можливих способів її вирішення розроблений 
в праці |71, інший - названий ,гранична задача Алексідзе для рівняння Лапласа" - пропонуємо вашій увазі. 

Постановка задачі |8). Предметна модель задачі відновлення потенціалу за МГПСТ - проста модель Зем- 
лі як абсолютно твердого тіла, близького до тіла обертання, що рухається рівномірно вздовж орбіти, обертаю- 
чись навколо осі з постійною кутовою швидкістю (без прецесії і нутації). Якщо у  - обмежена область про- 
(з) 


23 д - - » і б . 
стору К//, зайнята масами Землі, у - необмежене доповнення до У , вільне від тяжіючих мас, бу - фізи- 


чна поверхня Землі - границя множин Уу |і у" , то в прямокутній декартовій системі координат Ох х,хз з 
початком у центрі Землі, осі Ох, з Ох, якої лежать в екваторіальній площині, а вісь Ох; співпадає з віссю 
обертання, потенціал сили тяжіння мас всередині Землі м(2) С є у згустиною ам (2)- с(ає такий: 
с(С О(х) хе у 
УРНИ КА ЗУРИ (9 їн (1) 
жено 0,  хєу 
кі 


де / - гравітаційна стала, О(х)- 0.5 (х, б хо - потенціал центробіжної сили, (2 - модуль вектора ку- 
тової швидкості Землі. Напруженість поля (значення МГПСТ за (41) дорівнює 


Коорутоу зт о. г 


«ЛОдх, п 





дхі (х) -- т са 
де --Лоо- -сов(п, Х, , К -1,2,3 - напрямні косинуси одиничного вектора п(х) внутрішньої нормалі до 


дп 


еквіпотенціальної поверхні ПУ (х): У(у)- СХ, яка проходить через точку Х (рис. 1). 
Нелінійна гранична задача Алексідзе для рівняння Лапласа: необхідно знайти функцію Ух), хє зе 
б з хо аа б те - б 
яка задовільняє / всередині необмеженої замкнутої області у -у обу рівнянню Лапласа 


ДУ) -0, хє у", якщо в будь-якій точці ляпуновської границі бу області і в нескінченно віддаленій точці 


2 
| дУУ(х 
вона задовільняє умовам: УЗ що зах), хєдбу, УУ(х)-» 0 при м -»00, де «(х) - задана 
кеі Хі 
неперервна функція. 
Гармонічну в області у" функцію Ух), хе у" природньо шукати як потенціал простого шару |9) типу 
(1) (без урахування центробіжної складової О(х)) з невідомою густиною с(х), хе бу (інтегрованою, хоч 
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загалом може бути більш гладкою), поширеною на поверхні Ляпунова бу. Рівняння, з якого відновлюють не- 


відому густину сх) за заданими на поверхні бу значеннями хх) МГПСТ, виведено, виходячи із зо- 
браження: 





8 "(х) о що со5(р, 4)45, 45. (2 
Іл? з 3, хо і п 
де одиничні вектори р і 4, спрямовані відповідно з точки Х в точки є і 7, що пробігають по поверхні бу, 


хо- 


і Я: 


, 4 я с0б(х,4)- ет) 





мають вигляд р; - со5(р,х,)- , а кут між самими векторами - 
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Аналітичні властивості функції МГПСТ б характеризує таке твердження (10): МГПСТ не задовільняє 
рівнянню Лапласа в жодній точці області у" , що випливає з подання (2). 

Крім того, функції МГПСТ «(х) простого шару властива неперервність, що випливає з наступної леми. 

Лема 1. Модуль градієнта потенціалу простого шару є неперервною функцією точки Х є бу , яка рухаєть- 


ся по поверхні бу Ляпунова. 
Лема 2. МГПСТ простого шару, поширеного на сфері радіуса р з одиничною поверхневою густиною 


с(х) -1, х є ду, визначає вираз 
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Це означає, що функція ах) розривна, і має розрив неперервності при переході точки Х через повер- 
со8(п, 1 
хню бу. Доведено, що величина цього розриву при с(х )- 1, хе бу за умови сон за 2 


» реоеЇ 
і 2 2 
для сфери дорівнює 8; (х)- 26 (х) -3/4, х є бу , що узгоджується з лемою 2. 
Розв'язання задачі. Якби на поверхні Землі (рівняння якої є заданим), крім значень МГПСТ (х), хєду 
ії внутрішньої нормалі т(х) до бу, вимірювали напрямок градієнта п (х), задача відновлення потенціалу 


Н - З з 7, Й 5 
сили тяжіння У(х), х є у" звелась би до розв'язання зовнішньої задачі Неймана для рівняння Лапласа 


ТО 1 хе ду, ТОД 





УЗУУ(х) - 0, хє у", -»00, (3) 








3 3 
зауважуючи, що Ф(х)- о - б - я(х)соє (п, т)- ' Ус, сов(п,т)- 5 сов(п, х, )со8(х, ,т). 
т т Езі Каі 


Але напрямок нормалі л(х) (рис. 1) невідомий через виняткову складність і вартість вимірювань. Граничні 


дані задачі відновлення потенціалу за значеннями МГПСТ «Сх)- - Уа), у рамках прийнятої моделі 


описує вираз 
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або |6) г С) в? (х уз со5( п со) -я ооо , через те, що 200) 
ке дху й Хз 


-0. 





Введімо нормальний потенціал, що генерується фіктивними масами, що дорівнюють масам в Уу , але роз- 
ташовані в певному сенсі "нормально" в деякій іншій області у, простої геометрії з границею дуб» яка не ду- 


же відхиляється від земної поверхні бу . Подаймо потенціал притяжіння у вигляді суми Ух) ш С(х)ч т(х) 
З 


нормального С (х) і збурювального т(х) потенціалів, завдяки чому збурювальний потенціал описує відхи- 
лення реального розподілу масв у від нормального. Нехай У(х) - одинична внутрішня нормаль до повер- 
хні 00, Су) «С, а у(х) - Уц(х) - модуль градієнта нормального потенціалу, і задані напрямні ко- 


синуси собу, Хі ), соз(х, ,т) внутрішніх нормалей У(х), т(х) до поверхонь ду, і бу , аразом з нимиї 


3 
со5(у, т)а У соз(у, Хі )соз(х, ,т). За таких припущень відновити потенціал притяжіння УУ(х), хе у 
ко 
можна з граничної задачі (3) шляхом обчислення послідовних наближень УЗ), Кк « 0.1,2,...,00. Алго- 
ритм такий: 
1. за знайденими з попереднього і-го кроку наближеннями сок(п,, х, ), К с 1,2,3 напрямних косинусів 


нормалі п(х) визначаємо на границі бу за формулою (4) і -- І-ше наближення сили тяжіння (7) 
У 
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з до) я 


зе в(х 2. соз(п, ху )- 8і "(чле дк 





соб(п,, т) У сов(п, ,х, )соб(х,, т), хєду, 
Ка 
Фа (2) В він (хсоз(п, , т) -б у(х)сов(у, т), х єдбу; 
2. знаходимо розв'язок зовнішньої задачі Неймана для рівняння Лапласа 


па з «Фах ) хеду, По ) 


(х - х є бу, розподілених з на поверхні ДУ: 


АТ а (х - 0, хеу», 


як потенціал простого шару мас неперервної густини д.,, 





стрясь, (5) 





Т(х)- еПези За хау": 6) 





3. невідому густину ше з нелінійного інтегрального зе Фредгольма 2-го роду |6Ї: 


а Ук х б Бл (2345, 2 2Фуц (2), хе ду, т? 





1 д 1 1 со5(и,т з РО 
де К(х,б)- - ) сок(и, т)- 5 сок(т, х,)ці-яі це уижевхоі. 
2я дт,|х -4|  ш2я р-й «о ре 
4. Розв'язавши рівняння (7), наближено обчислимо з використанням (6) похідні потенціалу притягання 
| оУЗЗ() дц(х) дТ. (х) ОТ га 
У) ре ( ) се ( и з ) ЕЕ б, і Е- ба (2)45.., 1 лоб -1, 2. 3; (8) 
Ох, Ох, Ох, дх, Ал х-є | 
З (8) очевидно, що похідні наближень збурювального потенціалу визначаються у внутрішніх точках області 








я -- . . - . . . г. 
у", анаїї границі ду значення похідних з (8) неможливо знайти через неінтегровні особливості у підінтег- 
ральних функціях. Для продовження обчислень слід знати значення похідних збурювального потенціалу саме 
на границі бу , і для їх обчислення слід передбачити спеціальну регуляризацію інтегралів (8). 


Якщо цю операцію виконано і знайдено значення похідних (8) в точках Хх є бу , обчислюємо наступні наб- 


| РИНО ВЕ 
соло ) - б ТЕЗ АНЕ хє бу, рн (х) ж з(х ум ,Х :)- зх у-2- аа ,хеду 


кі дх 


лиження (крок 1 алгоритму) 


Фо (0) З вно (0)со5(лу т) У(х)соз(у, т) х є ду, 
далі знову розв'язуємо граничну задачу (5-7), наростивши індекс і -- 1-го наближення збурювального потен- 


ціалу Т.,, (х) і наближення У) - (Сх ут 


а ) потенціалу притяжіння (кроки 2-4) і т.д. 


Обгрунтування задачі. Заміна коректної задачі (5) розв'язком граничної задачі (6)-(7) і збіжності набли- 
жень УЗ) до потенціалу притяжіння УУ(х), хе ді строго обгрунтувана на прикладі близької задачі |4). 
Однозначчність її розв'язку доводиться теоремою єдиності розв'язання задачі Неймана для рівняння Лапласа 


за МГПСТ через потенціал простого шару, зведеною до доведення збіжності чисельних наближень Ух), 

К « 0,1,2,...,00 функції Ух), хе зо Її ж легко довести, виявивши збіжність послідовності т, (ж): зі збі- 

жності Шт Т, (х) -» т(х) випливає збіжність Ух) УУ(х), х є у", і збіжність до своїх границь будь- 
К-зоо 

яких інших наближень, що однозначно визначаються за Т, (х). Аналогічна стратегія справедлива і для даної 

задачі. Зокрема, справедливі наступні теореми |9). 

2 
ут(х) 
дО 
УС(х) 


розв'язків т, (х)) граничних задач (5) збігається до збурювального потенціалу т(х) області у . 


Теорема 1. Якщо величиною є (х) можна знехтувати порівняно з «(Х), то послідовність 


За умови Їсох(у, т)| щ- І (напрямок внутрішньої нормалі у(х) до еквіпотенціальної поверхні ди, збі- 


гається (або протилежний) з напрямком нормалі т(х) до поверхні Землі ду) задача зводиться до зовніш- 
ньої задачі Неймана для рівняння Лапласа. 

Формули (8), справедливі для внутрішніх точок області у , можна поширити на граничні точки Х Є бу ,на 
що вказує наступна теорема. 

Теорема 2. За неперервної на границі бу функції густини потенціалу простого шару (6) граничні значення 
частинних похідних потенціалу 1-го порядку дорівнюють 


от (х) | 1 АГ о по 1 й 
б і Ї ян 5(2)45, - 5 б) сок(х, і; 


Ця формула непридатна для практичного обчислення похідних (через складність земного рельєфу), тому 
замість неї варто використовувати еквівалентну формулу 


ПО бал, ко б(бвовоху т, ). 


Це нелінійне рівняння, як і попередні, з теорем 1 і 2; за фіксації геометрії контактної поверхні дУ(), зада- 








2 ; А - : го Й 
НОЇ на класі Ляпунова С! 6 ) задача знаходження густини потенціалу простого шару стає лінійною І має 


однозначне розв'язання. 
Щодо невідомої густини потенціалу простого шару виведено нелінійне інтегральне рівняння сили тяжіння: 


1, ; 1 2 
47 (хз роми с(Є 345, чн пре ГП о а соз( р, 4)45, 45. 2: 8 (х), х є ду. (9) 


Його розв'язок с), хе бу еквівалентний розв'язку задачі Алексідзе з граничними даними на поверхні 











ду ре- | дуду 


- б . Ж . 
ду Ляпунова, оскільки за будь-якого вибору густини потенціал простого шару задовольняє в області у рів- 


нянню Лапласа, а знайдене з (9) значення густини забезпечує виконання граничної умови. Питання розв'язно- 
сті задачі Алексідзе редукується до вияснення умов існування, єдності і стійкості розв'язку рівняння (9). 
Це рівняння можна спростити до такого вигляду: 


о'(х)ч би П б т со5(р, 4)45, 45; б 28 (х), хєду. (10) 


дусу|/ - -1 


Розв'язки рівнянь (9, 10) допомагають визначати не лише потенціал Ух), хє з а й значення МГПСТ 





у й з б - . б Е Р 
в будь-якій точці необмеженої області у . Останні можна обчислити як за виразом (27), так і за зручнішою 


для обчислень формулою 
2 


з"(х)е ново, | (1) 


яка не вимагає обчислення двократного інтеграла. 

Чисельні алгоритми. Зведення задачі Алексідзе з граничними даними на поверхні Ляпунова до розв'я- 
зання нелінійного інтегрального рівняння сили тяжіння (9) чи (10) дозволяє легко вивчити питання розв'язності 
та єдиності її розв'язків, і ефективно знаходити чисельні наближені розв'язки у випадку областей складної фо- 
рми. Дослідження питання розв'язності задачі Алексідзе зводиться до вияснення умов коректності рівняння 
(9). Конструктивні умови єдиності, існування та стійкості розв'язку задачі, подані у вигляді відповідних теорем, 
дозволяють вказати такі ітераційні схеми для її обчислення. 

Алгоритм 1. Для визначення густини потенціалу простого шару придатний ітераційний процес 








Со (2) - 2(х) ори (х) -0,хебу, бу (2) -Ої, (х)ч о (2) кору (2) ЗА Їду сої па Осо , 





1 , 
де А, сь С - ост, )- Гкуь ба, у», (235, ; К (о босі)з 2 ст) сові 45, ; 
бу 4л ду о (2) Їх лю п 


Мосау)я б (оР09- Ва). 


с (х) 


Алгоритм 2. Визначенню послідовних наближень густини слугує ітераційний процес 


10 (2 з х(х), еру (9 -0, хєду, Фін (2) А, (х)ч о (2) бу (2) -дА, Їсть, (9 пе Осо , 
ле Ас», 9 - ь (ват, )- Гк, (зба, у», (235. , ь (зоу)- (27(9)- Б (сет )- В оо» )/ ау (х) 


Алгоритм 3. Обчислення послідовних приростів густини з лінійних інтегральних рівнянь 2-го роду 


сна Х)Я Гк. (о балон (Є 5, бе ь(» ві 
ду 


щодо приростів С», (х) зобі (х)- ор (х), паб; Со (х) ш «(х.), ори (х) -0, х є ду густини. 
Щоб замкнути теорію вирішення задачі Алексідзе, лишається дослідити питання скінченновимірної апрок- 


Висновок. Сформульована нова нелінійна гранична задача Алексідзе для відновлення потенціалу за зна- 
ченнями МГПСТ і вказано алгоритм її розв'язання як потенціалу простого шару (6), густину якого практично 
відшукують з рівняння (7). Знаходити густину з еквівалентного нелінійного інтегрального рівняння сили тяжін- 
ня (9) недоцільно, але воно зручне для вивчення умов коректності її постановки з граничними даними на класі 


(2) (у 


Ляпунова С 4 Задача Алексідзе з граничними даними на поверхні Ляпунова редукована до розв'язан- 


ня двох еквівалентних нелінійних інтегральних рівнянь (9-10), що описують функцію сили тяжіння, У цих реду- 
кціях вона коректна у на парі банахових просторів, до яких належать вхідні дані і шуканий розв'язок. 
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